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Abstract. Dynamics ef neutral modes for fractional quantum Hall states is investigated 
for a quantum point contact geometry in the weak- backscattering regime. The effective 
field theory introduced by Fradkin-Lopez for edge states in the Jain sequence is generalized 
to the case ef propagating neutral modes. The dominant tunnelling processes are identified 
also in the presence of non-universal phenomena induced by interactions. The crossover 
regime in the backscattering current between tunnelling of single-quasiparticles and of 
agglomerates of p-quasiparticles is analysed. We demonstrate that higher order cumulants 
of the backscattering current fluctuations are a unique resource to study quantitatively the 
competition between different carrier charges. We find that propagating neutral modes are 
a necessary ingredient in order to explain this crossover phenomena. 



PACS numbers: 73.43. Jn, 71.10.Pm, 73.50.Td 
1. Introduction 

The Fractional Quantum Hall Effect (FQHE) represents one of the most important examples 
of strongly correlated electron systems [Hill [3]. Due to the purely two-dimensional electron 
motion a wider class of symmetries, for the exchange of identical particles, are admissible in 
comparison with the three-dimensional case [Ij. This implies the existence of quantum fluids 
with unique properties, such as fractional charges and fractional statistics, and eventually 
different internal topological orders [5] at a given filling factor [6]. In order to discriminate 
among these different theories it is useful to investigate and characterize the edge state 
dynamics [7j. Indeed, one may achieve unique information about the bulk properties via the 
holographic principle [8], by studying the low energy gapless excitations at the boundaries 
of the Hall bar [9l[l0l. 

Edge states are described in terms of bosonic modes of a Chiral Luttinger Liquids (xLL) 
For the description of the Laughlin sequence [1] with filling factor v = l/{2m + 1) (m G N), 
only a single charged mode is necessary. Here, theoretical predictions of transport properties, 
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as current and noise, through a Quantum Point Contact (QPC) were proposed fTH and 
confirmed by experimental measurements [Hi [15]. It was demonstrated that, in the weak- 
backscattering regime, the tunneUing carriers have a fractional charge e* = ve (e electron 
charge), giving the first confirmation of charge fractionalization [1]. Further experiments [H] 
attempted to demonstrate also the anyonic nature of these excitations |17j . 
For the Jain's sequence [18] with ^ = p/ {2mp + 1) and p e N one needs additional 
bosonic neutral modes in order to fulfil the hierarchical structure [19]. As suggested by 
Wen and Zee [20], one can introduce p — 1 neutral modes in addition to the charge one. 
The model predicts the expected value e* = ve/p for the charge of the fundamental quasi- 
particle (qp) excitation, however, due to the increasing number of fields, for increasing 
hierarchy p, the conserved currents grows indefinitely [2T1 [22], creating some problem 
of consistency [23]. Some predictions of this theory have been experimentally verified 
by investigating the transport properties in cleaved-edge overgrown devices but, at the 
same time, experiments shed light on the limits of the model [2H [25] . In particular 
the behaviour of the tunnelling current suggests that the dynamics of neutral modes is 
less relevant with respect to the prediction [12], [26] . These problems triggered additional 
theoretical investigations devoted to clarify the reasons of the discrepancies [27]. Among 
them, it is important to mention the model proposed by Fradkin-Lopez (FL model) [23], 
which introduces only two neutral topological modes (one for infinite edge length |28j). 
instead of the above dynamical p — 1 modes. These modes have zero velocity, then 
they do not contribute to the dynamical response, having only a role in enforcing the 
appropriate statistical properties. This model predicts also the correct value of the single- 
qp charge in agreement with the first experimental observations on tunnelling through a 
QPC [29]. However, successive experiments [30l [31] have shown unexpected behaviour in 
the conductance and, more intriguingly, a value ~ pe* for the carrier charge at extremely 
low temperatures, larger then the single-qp. These observations are not explainable within 
the framework of the original FL model. Recently we suggested that these anomalies can be 
understood by generalizing the FL model assuming propagating neutral modes [32]. This 
explain the peculiar behaviour of the experiments and the crossover regime between two 
different carrier contributions: the single-qp with charge e* and the p-agglomerate of qp 
with charge pe* . 

The role of neutral mode dynamics has been also addressed for other FQHE states [33 l [M l [35] 
with emphasis on the concept of finite velocity propagation [35l [36l [37] . It is believed that 
Coulomb interactions strongly increase the charge mode velocity in comparison to the neu- 
tral one such that the two modes energy cut-off bandwidths differ consistently [35], [36] . 
Recently, several theoretical proposals on the possibility to detect the neutral mode velocity 
have been formulated. They range from the thermal Hall conductance [2T1 [22] and the res- 
onances in tunnelling through an extended QPC [38], to the analysis of Coulomb blockade 
peaks [39], [40], [41] . Unfortunately, from experimental side the evidence of neutral modes is 
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often missing or weak and quite indirect. 

The work of the present paper is motivated by these open questions, with the main task of 
giving evidence of the dynamics of propagating neutral modes in transport properties. In 
order to produce results that can be tested experimentally we choose the standard geometry 
of a QPC for the Jain series accessible to nowadays experimental capabilities. A possible 
renormalization of the dynamical response, due to external interactions [l2lll3llllj, is also 
taken in account, we will show that this introduces non- universal behaviours. 
In addition to current, we will analyze also current fluctuations [15] in order to properly 
describe tunnelling of different charges. The available theory of current cumulants for a 
general tunnelling Hamiltonian [16], combined with the recent successes in measurements 
of current cumulants in QPC [UlllH] and other structures [l9l[50], gives a firm background 
and a strong motivation for this analysis. Current cumulants |5ll [521 [53] contain indeed 
several information on the carrier charges, and may discriminate between models based 
either on a single carrier or on many carriers, considering a direct comparison with the 
available experimental data on QPC in Hall bars [301 [HH [53, [56] . 

We first find that neutral mode dynamics may be detected in the current behaviour of 
single-qp for low enough temperatures. We also demonstrate that weak interactions are 
more favourable in order to detect the effects. 

We point out that the most relevant effect, signal of the presence of neutral modes, is the 
possibility to have different tunnelling excitations by varying temperatures and voltages. 
We prove the existence of a crossover regime which separates the tunnelling processes of 
p-agglomerates, at energies lower than the neutral modes bandwidth, from the standard 
tunnelling of single-qps at higher energies. We demonstrate that it is possible to find voltages 
and temperatures regimes where the p-agglomerates are detectable, although in presence of 
interactions and larger temperatures these effects cannot be easily visible. For this reason 
we analyse noise and higher cumulants of the current fluctuations, whose behaviours can 
give clear information on the crossover regime. 

In the shot noise regime the Fano factor is a proper quantity to detect the effective charge 
of the dominant carrier involved in tunnelling. We observe, as a function of the external 
voltage, a double plateaux structure clearly indicating the contribution of both the p- 
agglomerates and the single-qp excitations. In order to avoid the thermal regime, present 
in noise measurements, we consider higher order odd cumulants, such as the skewness. 
These, indeed, are less sensible to temperatures and show plateaux structures that can be 
interpreted in terms of dominant contributions of different tunnelling charges. We show that 
the comparison between different current cumulants behaviour can be a unique resource to 
correctly interpret the physics of this crossover phenomena [301 ESI EEl [57] . 
The outline of the paper is the following. In section [2] the generalization of the FL 
model to the case of finite bandwidth neutral modes is presented and the field operator 
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construction of the edge excitations is derived with particular care of the monodromy 
condition. The scahng dimension is then discussed in order to find the dominant tunnelhng 
processes through the QPC In section [3] we discuss analytical formulas for transport 
properties via Keldysh technique at lowest order in tunnelling. The relation between current 
cumulants and backscattering current is presented. Results are discussed in Sec. |4]both on 
the backscattering current and on noise and higher current cumulants. Conclusions are 



presented in Sec. |5l The [Appendix A| contains the detailed derivation of the field operator 



construction of the edge excitations. 
2. Model 

We start by recalling the description of the edge states dynamics at the boundaries of an 
Hall bar in the Jain sequence [18] with filling factor u = p/{2p + 1) (p G N). Within the 
minimal FL model [231 EHl [32] , edge states consist of two counter-propagating bosonic fields: 
a charged mode ip'^{x) and a neutral one <f^{x). The corresponding Euclidean action for a 
given edges side is [23], [32] (from now on h= 1) 

1 f f^^"^ 











r.L/2 
'-L/2 
■.L/2 
-L/2 



+ J J d^^'^ix, t) {-idr + Vnd^) ip''{x, t) , (l) 



where /3 = {kBT)~^ is the inverse temperature, and L is the edge length. We will consider 
the case L — oo neglecting correction induced by the finite size. Following [32] the modes 
propagate with different, but finite velocities: Vn for the neutral mode and for the charged 
one. The velocity of the charge mode can be strongly increased by Coulomb interactions so 
we can safely assume f„ <C Vc [35] . 

As we have already shown in [32], neutral modes with finite velocity affect qualitatively the 
transport dynamics and could explain the observations of recent experiments [30]. We will 
investigate this general case focusing on the main effects induced on transport properties. 
Notice that neutral modes with finite velocity differ from those described in the original 
FL model, where they are topological (f„ = 0), therefore we denote the present model as 
Generalized Fradkin-Lopez model. 

As clearly emerge from the action ([T]) the two modes at the same edge propagate in opposite 
directions. Notice that, differently from what happens in the hierarchical model [20], where 
one needs p fields to describe each edge, in this minimal model the number of fields is always 
two, independently on the value of filling factor of the Jain series [23l [28] . 
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The bosonic fields in ([T]) can be represented in wave-number components as 



fc>0 

where k = 27is/L (s G N) are the wave vectors, I = c,n denotes the charged and neutral 
modes and i>c = z/, i/„ = 1. The parameter represents the ultraviolet wave-number cut- 
off that, together with the mode velocities, determines the energy bandwidths = Vc/a, 
= Vn/a- In the following, ojc will be considered as the greatest energy scale with 
ujn <^ uJc- The bosonic creation and annihilation operators obey standard commutation 
relations [a\k) , a''' {k'^] = 6i^ii6kk'- They ensure that the fields satisfy [20], [23] (for L ^ oo, 

[V9^(x), ip^iy)] = i7ri/sgn(x - y) ; (f'^iy)] = -Z7rsgn(x - y). (3) 

The electron number density along the edge is 

p{x) = ^d.y^^x) (4) 
and depends only on the charge mode. 

2.1. Quasiparticle operators 

To complete the description one needs an expression for the operators creating excitations at 
the edges. To be as much general as possible we will discuss both the single-qp excitations 
with charge e* = e (zz/g) (e the charge of the electron) and multiple-qp excitations with 
charge me* with m > 0||| We refer in the following to these as m-agglomerates [32] . 
There are three main constraints that the m-agglomerate operator '^^™-\x) has to 
satisfy [321 [60l |6T] 

(i) Charge: multiple excitations must have a charge of me*. This implies appropriate 
commutation relation with the electron density 

[pix),¥^\y)] =-m(^^^5(x-i/)vl>(-)(|/). (5) 

(ii) Statistics: excitations have to satisfy fractional statistics 

¥'^\x)¥'^\y) = ^M(y)^M(x)e-^^'"''*^''(^-^) , (6) 
with the statistical angle [TTl |23] 

VP P J 

The multiplicity of A; G Z takes into account the 27r periodicity of the phase. Note that 
this expression, given in the FL model, is in accordance with the standard hierarchical 
construction [20l. 



I The operators with m < can be easily obtained from the operators with m > by exploiting particle-hole 
conjugation. For m = we have neutral excitations that do not impact on the current properties. 



Neutral modes edge state dynamics through quantum point contacts 6 

(iii) Monodromy: the phase acquired by any excitations in a loop around an electron must 
be a multiple of 27r |62l 163] . 

Using the bosonization technique the operator associated to the m-agglomerate can be 
expressed as an exponential of linear combinations of the bosonic edge fields [321 [6T] 

\^("^)(^X) = — gi[amV''(a;)+/3m</3"(a;)] ^ ^g-^ 

\/2-Ka 

Here, JF^™) represent the so called Klein factors which will be discussed later, and the 
coefficients am and f3m are determined by imposing the fulfilment of the above three 
constraints. A detailed discussions of the possible form of these coefficients can be found 
in [61] for u = 2/5, and more generally in Appendix A for 1 < p < 6. Hereafter we simply 
quote the main results. The relation ([5]) fixes a^- Indeed, using the commutation rule ([3]) 
one finds directly 

m 

am = — ■ (9) 

p 

The second point (ii) is fulfilled by imposing the relation on the neutral par 



Pm=\lm^ l^l + i^ +2A; . (10) 

As one can see jSm still depends on the free parameter in ([7]), showing that, for a given 
m, there is a family of different excitations each with the same internal statistics. From the 
square root in (fTOj) it follows 



k ^ km,in, — lilt 



ni^ f 1 



(11) 



where Int [x] represents the integer part of x [Jj]. 

The last requirement (iii) reduces even more the possible values of k. As shown in Appendix 
A, for p < 6, it is possible to find the following form 

k = ^^^^{q' - s') + dip + l)(g - s) , (12) 
where g G Z, and s and d are integers that identify the m-agglomerate for a given p 

m = sp + d, with < d < p — 1, and s > . (13) 

Inserting the expressions ([T^ and f[T^ in flTU]) one obtains 



(14) 



§ The over all sign in the expression for (3m is arbitrary, the physical results on scaling dimension are not 
affected by this choice. 

II Note that we changed definition of the coefficients am and (3m and the sign of k with respect to [32]. 
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where, for sake of clarity, we explicit the dependence on q. We can now write the final 
expression of the field operator, for a given q, within the family of the m-agglomerate 

.l.(m,q)(^^^^^ _ •^^"''''^ ^i[(s+|)y-(x,t) + Vp(p+l)(q+|)y"(x,t)] _ ^^^^ 

We conclude by commenting on the role of the Klein factors JF^™'"?). They are "ladder 
operators" that change the particle numbers of the (m, g) -agglomerate type on the 
edge [m [65]. They ensure the appropriate statistical properties between excitations on 
the same m-family, with different q [61] or on different edges [66]. In the following we will 
omit them since, apart their straightforward ladder operation they do not play any role in 
the sequential tunnelling regime treated in this paper [66l [67] . 

2.2. Scaling dimension 

To investigate the relevance of different m-agglomerates it is useful to consider the behaviour 
of the two-point correlation function at imaginary time [3] 

^(-''?)(0,r) = (T,^("^''')(0,r)^t(-''')(o,o)), (16) 

where is the imaginary time ordering operator and (...) the thermal average with respect 
to the free action ([T]). To evaluate ( |T6l) one can use the bosonic representation of the 
excitation field f[T5l) ending up with to two distinct averages on charged and neutral modes 

27ra 

^ J_gam»c(0,T) /3^{g)©„(0,r) 

27ra ■ ^ ' 

Using standard identities for boson operators the above averages are represented as standard 
bosonic propagators (/ = c, n) 

A(0, r) = (TV(0, r)^'(0, 0)) - ((^'(0, 0))^) . (18) 

Since the Euclidean action ([1]) is quadratic, Vi^O, r) may be easily computed. For the present 
discussion we only need the T = limit. One has [3l [68] 

Vi{0,T) = -uiln{l + iJi\T\), (19) 

where it is evident the role of the bandwidths ui as cut-off energies. Substituting this 
expression in (fT7|) one obtains the zero-temperature correlation function 

1/1 \ ''"m / 1 \ f^mil) 

^(-.'?)(0, r) = — —— —— . (20) 

^' ^ 27ra \1 + UJc\t\ J \l+UJn\T\J ^ ' 

We can now compare the scaling properties between different (m, g) -agglomerate operators. 
First we recall that the relevance of a given excitation is obtained by considering the scaling 
dimension Am(Q') defined as half of the exponent of the correlation function fl2Ul) in the long 
time limit, at T = 0. From ([20l) one has, for r — > oo [31 [26] . 

e?(™''')(0,r) oc Irj-^^'-^''). (21) 
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The most relevant operator is then the one with the minimal scaling dimension, this 
corresponds to a dominance at long times and consequently at low energies. 
The scaling dimension of (m, g)-agglomerates is easily derived from fl20|) and fl2T]) 



u [ s + - ] + pip + 1) g + - 
Pj \ P 



(22) 



It depends both on charge and neutral modes contributions. 

We are now in the position to determine the most relevant operator at low energies. Let us 
start to consider Am{q) for a given m-family [s and d are fixed by m according to (fT3ll ). 
It is easy to see that the most relevant operator among different g's is given by g = for 
d/p < 1/2 and q = —1 for d/p > 1/2. Concerning the d values there are two possible 
classes: d = which corresponds to s multiple of j9-agglomerates (m = sp) and d > 1. For 
the first class [d = 0) the most relevant operator corresponds to s = 1 with m = p, denoted 
as p-agglomerate with 

Ar = r<l- (23) 

For the second class {d > 1) the minimum scaling is reached for d = 1 and s = 0, namely 
the single-qp with 



1 






1 








> - 


2 


.P \ 


pJ 


2 



(24) 



Comparing fl23l) and fl24|) we conclude that the p-agglomerate is always the most dominant 
operator at least for the discussed case with 1 < p < 6. This result covers the standard 
scaling region with energies much smaller than both the neutral and the charged bandwidth 

In view of future analysis of transport properties it is also important to discuss the 
intermediate energy region with u;„ <^ <^ lJc, despite due to the finite energy values 
the standard scaling argument has to be taken with more care. In the intermediate regime 
t^n ^ |t|^^ ^ ujc, the neutral modes in (!20!) are saturated and only charged modes contribute 
to the dynamics. It is than convenient to introduce an effective scaling dimension 

^" = ^ = 1^. (25) 

that depends only on the m-agglomerate charge. At these intermediate energies we then 
expect the dominance of the single-qp with m = 1, indeed A^''^ < A^^ for any p > 1. Then, 
the presence of a finite cut-off u;„ discriminates among two different regions with a crossover 
from energies (e.g. voltages or temperatures) smaller than a;„, where the p-agglomerate 
dominates, to energies larger than Un where the single-qp could become more relevant. 
Before to conclude we would like to discuss the robustness of the above results in the presence 
of possible interaction effects due to additional external degrees of freedom. It is well known 
the possibility to have renormalization of the dynamical exponents induced by a coupling 
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with external dissipative baths . Different mechanisms of renormahzation were proposed 
ranging from coupling with additional phonon modes [l2], interaction effects [7Ql|7T], edge 
reconstruction |l3lllll[72]. In order to qualitatively take in account these effects, we follow 
[l2] , considering baths given by one dimensional phonon-like modes, coupled to charged and 
neutral modes at each edges. As shown in detail in [32] this interaction renormalizes the 
imaginary time bosonic propagators (JT8l) as 



with gc and gn interaction parameters. Note that gc correspond exactly to the factor F 
defined in [12]. Here, we also assume renormahzation of the neutral modes described with 
the parameter gn- Within this model it is gn,c ^ 1 with (yf„ ^ = 1 for the unrenormalized 
case. Note that the above renormalizations do not affect the field statistical properties, 
which depend only on the equal-time commutation relations, i.e. the field algebra. In our 
discussion we simply consider the renormahzation factors gc and gn as free parameters. 
Following analogous calculations, as done before, one can easily shown that, at low energies, 
the p-agglomerate is still the most dominant operator under the condition 



otherwise the single-qp will dominate. Note that, from fl27|) . strong charge renormalizations 
may destroy the dominance of the p-agglomerates in favor of the single-qp. However, 
the necessary values of renormalizations must be so high, especially if some neutral mode 
renormahzation is also present, that this probably would not easily happen in real cases. 
This explain why we expect that the p-agglomerate should be the most probable tunnelling 
entity at low energies in experiments. On the other hand, at intermediate energies 
tUn -C <^ Uc, the relevance of the single-qp it is not affected by the presence of interactions 
and it is then more robust. 

2.3. Tunnelling Hamiltonian 

Since the first measurements of transport properties on Hall bars through QPC [T H flS l [29], 
tunnelling experiments revealed to be a precious tool to investigate edge dynamics and 
internal properties of excitations. In this paper we consider a single QPC with tunnelling 
between right and the left side edges. 

In the previous analysis we discussed only one edge, now we have to extend it to the case of 
two different edges. We can simply set the curvilinear abscissas on the edges such that right 
(j = R) and left (j = L) edge fields V5j(x) have exactly the same field algebra as presented 
before for the single edge. 

Tunnelling of an m-agglomerate is described by the standard Hamiltonian 




(26) 




(27) 



(28) 
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where '^^p\xo) indicates the annihilation operator of the m-agglomerate on the edge 
j = R,L at point xq. This operator is defined in f[T^ for a generic edge. The coefficients 
and (3m will be chosen to have the minimal scaling dimension (122|) which, as we discussed 
before, corresponds to the most relevant operator within a given m-family. Without loss 
of generality we fix the point of tunnelling at xq = 0. The factors tm are the tunnelling 
amplitudes, they depend on the detailed structure of the edges, not known in the framework 
of effective theories, and on the precise geometry of the QPC, difficult to be experimentally 
controlled. In the following we will define them in terms of the single-qp tunnelling matrix 
element ti = t such that \tm\ = /tm|t|, where the real factors Km, represent the relative 
weight of the m-agglomerate amplitude in comparison to the single-qp one (ki = 1). These 
parameters may be in principle fixed by fitting the experimental data. 
In the following sections we will calculate transport properties at lowest order in the 
tunnelling, assuming the amplitudes tm sufficiently small to fulfil the weak-backscattering 
condition for all possible m-agglomerates. In experiments one can appropriately tune the 
QPC with pinch-off voltage to always satisfy the above condition. Following the previous 
discussion on the most relevant operators in the different energy regime we will consider 
as tunnelling processes only the two dominant contributions: the single-qp and the p- 
agglomerate respectively. The total tunnelling Hamiltonian is then 

Ht=Y1 4"^^ = E '^PmP\0) + h-c. . (29) 

m=l,p m=l,p 

As already discussed in [32] this seems to be enough in order to explain the recent 
experiments in the Jain series through QPC [30] . 



3. Backscattering current and noise in QPC 

3.1. General Relations 

In the following, we derive general expressions for the tunnelling current and noise properties 
of the generic m-agglomerate. We will make use of the Keldysh formalism [671 ESI [71] without 
entering into the details and referring to the literature for further discussions. At lowest 
order in the amplitude, the possible tunnelling processes described by ( l29l) are independent, 
and can be evaluated separately. The backscattering current operator for an m-agglomerate 
is 

= tme* \tm e^'"'*^* ¥^'>\t)¥^\t) - h.c.l , (30) 



where we adopted the interaction picture with the tunnelling Hamiltonian (12811 as interaction 
term. Here, for convenience, we introduced the effect of the bias V, applied to the QPC 
between the two edges, in the amplitudes tm by the appropriate gauge transformation 



+ ^ + pime*Vt m 
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The averaged backscattering current (^(^I^"^^'j'^^{t) can be written in the Keldysh formahsm 

as 

((/<'»)))^(i) = 1 Y^/Tk r/(™)((^)e-i/„-*'<'«')l \ . (31) 



Here, the index r] = +, — label the times on the Keldysh contour with r] = + for the forward 
branch {—oo,t) and rj = — for the backward one (t, — oo), the exponential term contains 
the integral of the tunnelling Hamiltonian in the interaction representation on the Keldysh 
contour J"^ dt. The average is with respect to the thermal density matrix of the unperturbed 
Hamiltonian at t = —oo and Tfi'[...] is the Keldysh time-ordering operator. 
Expanding the exponential term one obtains the perturbative series of the current in terms of 
the tunnelling amplitude. Inserting the bosonized form (fT5|) of the m-agglomerate operator 
and performing standard thermal averages as discussed in (JT71) . one expresses the steady 
(t —>■ oo) backscattering current at lowest order |67] 

//Iim)\\ ^ _ ^^f\ffm rZ'sHme*Vt') (32) 

where, Ifl'^ (t) are the Keldysh Green's functions relative to charged / = c and to neutral 
/ = n modes. They are related to the standard greater Green's function T>f{t) = 
(^'(t)(p'(0)) - (((^'(0))^) in the following way 

(Vf'iriltl) for 77 = 7?' 
1 Vf" iji't) for r] = —r( . 

Note that the right and left edges give an equal form to fl33l) . This is taken into account by 
the factor 2 in the exponents of ( l32i) . We recall that the greater Green's function are directly 
related to the temperature correlators introduced in (fTSjl applying an analytic continuation 
to real time [Ml EH] . 

Let us now consider the backscattering current noise correlator defined as 

((/(™)))2(t,t') = ^ {6I^"'\t)6l'^'^\t') + 6I^"'\t')6I^"'\t)) , (34) 

with 5I^"'\t) = I^^'^^t) — ((J^™)))^ (t) the current fluctuation. In the Keldysh framework 
one has [67] 

2 



(35) 

The second order in t can be now evaluated using the same procedure described for the 
current. In the stationary limit the current noise correlator depends on time's difference 
only, and the corresponding Fourier transform at zero frequency is [6j 
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(/(™)))^ = me*|((/(")))Jcoth 



(38) 



By exploiting the analitycal properties of the bosonic Green's function 0331) in the complex- 
time plane [69] 

Vr\t + .r^/3/2) = [Vr\t - .r//5/2)] * , (37) 

one can rewrite the time integral contour in ( l32i) and ( l36l) obtaining a relation, valid at the 
lowest order [671 [68], between the steady current (132|) and the zero frequency noise (1361) 

' Pme*V' 
_ 2 _ 

It is convenient now to discuss the connection with the standard perturbative approach of 
transport properties in terms of tunnelling rates r^'"^ The latter describe the probability 
per unit time to have tunnelling of an m-agglomerate between the two edges and can be 
easily calculated within Fermi golden rule at lowest order in t [32l [69] 

14-12^2 r+oo 

rM(E^) = / dt'e-*^™*'e2"^^=^(*') e^^^^"^*') , (39) 

where = me*V. As already pointed out by Levitov and Reznikov in [46] the tunnelling 
transport is described by a bidirectional Poissonian process. Therefore, there are two 
rates: the forward T^"^\Em) in (139|) and the backward T^'^\—Em), that corresponds 
to opposite tunnelling processes. They are related via the detailed balance condition 
T^"^\~E) = e~^^r^'")(i?) as imposed by the symmetry (137]) . From the knowledge of the 
rates characterising the bidirectional Poissonian process, one can directly calculate all zero- 
frequency current cumulants 



( (me*f (1 - e-^™^*^) T^'^^EJ k odd 

\(me*)'=(l + e-^"^^*^)r(-)(E^) even. ^ ' 

Note that k = 1 corresponds to the current fl32|) . while k = 2 gives the noise fl36|l . We can 
now discuss the linear regime f]me*V ^ 1, the linear conductance is given by 



G-M = hm ^\ = {me*f (3 r('")(0) 



,.Z-„ ,v <^^' 

and as expected from the fluctuation-dissipation theorem the thermal noise is ((/*■'"•' ~ 
2kBTG'-"'\ In the opposite regime (3me*V ^ 1 the noise has the shot behaviour [76j with 
((/M))^ ^me*((/('")))^. 

Relations similar to fl551) can be obtained between the zero-frequency current cumulants and 
the stationary current. Here, we simply quote the final result 

fcodd 



k 



[me 



[me*)^ ^ 



(/(™)))Jcoth 



l3me*V 



(42) 

k even. 



2 

Notice that all the previous relations are also applicable in the presence of renormalization 
induced by coupling with external baths, since they are essentially based on the general 
analytic properties (|371) of the thermal Green's functions. 
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Until now, we presented the cumulants for tunnelling current of m-agglomerates (see (125]) ). 
As already pointed out at lowest order, tunnelling of different agglomerates are independent, 
and can be evaluated separately. Considering the tunnelling Hamiltonian (l29l) with the two 
most dominant terms m = l,p the current cumulants will be the sum of the contributions 
of the two processes 

{{Itot)h= E^^^^'"^))'^- (43) 

m=l,p 

This independence will not hold at higher order in the tunnelling. 



4. Results 



4-1- Single quasiparticle current 



In this part we will analyse in detail the tunnelling current of single-qps with particular 
attention to the role of neutral modes. We remind that the single-qp operator is identified 
hj m = 1, d = 1 and g = in ([9]) and (|T4l) such that 



1 

1 + - 
p 



(44) 



The first ingredient necessary for the evaluation of the backscattering current (132!) is the 
rate (15^ that depends on the real time finite temperature Green's functions of the bosonic 
modes T>^{t). These functions are well known in literature |69J and can be evaluated by 
analytic continuation from zero-temperature imaginary time Green's function in fl26|) . Here, 
we simply quote the result also in the presence of interactions gn, gc 

\T{l + T/LJi-zTt)\^ 



V>{t)= giuiln 



r2 (1 + T/ui) (1 - ZUit) 



(45) 



with T{x) the Euler gamma function (from now on ks = 1). The evaluation of the rate will 
be done by performing a numerical time integration of ([31]). At zero temperature the rate 
can be also calculated analytically and the current is given by 



e* t 



27ra2 T{2i^gcai + 2gnl3t)Ei 
X iFi 



El 



2vgcCi\ 



El 



2gn(3l 2vg,ai + 2gn(3l 



(46) 



El 

with El = e*V and iFi[a,b,z] the Kummer confluent hypergeometric function. Note that 
in case of neutral topological modes, where a;^ = 0, the bosonic neutral mode propagator 
of (j^5|) becomes = 0. The corresponding current at T = is 



27ra2 T{2ugcal)Ei 



El 

UJr. 



(47) 
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Figure 1. Single-qp backscattering current as a function of tlie external voltage e*V/ujo, 
with ujQ =^ lO^^LUc, in bilogarithmic scale for ly = 2/5 and topological neutral modes 
i^n = 0. (a) Unrenormalized case gc — I', (b) renormalized case gc = 5.5. Different 
line styles correspond to temperatures: T/uiq — (solid black), 0.025 (dot-dashed blue), 
0.1 (dashed green), 2 (short dashed orange). The unit of the current is It — e|tp/(47r^ a^Wc). 
The zero temperature power-laws of the figures are given by V~^^^ and V^"^^^^ for the 
unrenormalized and the renormalized case respectively. 



which corresponds to the hmit c<j„ ^ of (146|) . showing the standard contribution of a single 
mode [691177]. 

To discuss the effects of finite bandwidth, we compare the current obtained for topological 
neutral modes in figure [1], with the one given in the presence of propagating neutral modes 
in figure O We will analyse in detail the case of filling factor u = 2/5, with af = 1/4 and 
(3f = 3/2. However, the results are qualitatively valid also for the other values of the Jain' 
sequence. We also consider the presence of renormalization. 

Figure [1] shows the voltage dependence of the backscattering current for topological neutral 
mode {un = 0) in the unrenormalized case gc = 1 panel (a) and in the renormalized case 
gc = 5.5 panel (b). At zero temperature, black lines, one recovers from ( 1471) the power-law 
behaviour 

(^(^lihtop)^^^^y2ua^al-l^ (48) 

driven by the charge exponent coefficient ai and with a charge renormalization gc- For 
the unrenormalized case figure IH^a) and for finite temperatures (coloured lines) the current 
exhibits a maximum around e*V ^ T with a linear ohmic behaviour for smaller voltages 
e*V -C T, and a power-law (HHl) at higher bias e*V > T. 
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Figure 2. Single-qp backscattering current as a function of tlie external voltage e*V/LOQ^ 
with LUo — lO^'^Wc, in bilogarithmic scale for v = 2/5 and propagating neutral modes with 
LOn = IQ^^ujc- (a) Unrenormalized case gc = g-n — 1; (b) renormalized case gc = 5.5 
and gn — 2. Different line styles correspond to temperatures: T/loq = (solid black), 
0.025 (dot-dashed blue), 0.1 (dashed green), 2 (short dashed orange), 20 (dotted red). At 
zero-temperature the power-laws scaling are: yf^/s+SS"-! for e*V <^ ujn and V^''/^~'^ for 
e*V 3> LOn- The unit for the current is as in figure [T] In the inset the linear conductance 
in unit of Gt = e'^\t\'^ / {A-k^o^loI) is represented in a bilogarithmic scale as a function of 
temperature. 



For sufficiently strong charge renormalization > /2v [g^. > 5 for z/ = 2/5), as shown in 
[D^b), the power-law exponent in (HH!) is positive and the current is always growing. Then, for 
finite temperatures the maximum is substituted by a smooth change of power-laws between 
the ohmic regime and the zero temperature behaviour fHHj) . 

We can now compare the above behaviour with propagating neutral modes as shown in 
figures [2](a) and[2](b). Let us start with the zero temperature case, here the current displays 
two different scaling 



oc 



V^^f^"!-^ e*V > iOn. (50) 



Neutral modes contribute to the dynamics only for voltages smaller than the bandwidth 
Un- The power-law exponent (H9ll is always positive because of (jHj) and gc,gn > 1- Only 
at higher voltages e*V ^ as in fl50l) . one recovers the standard FL behaviour of the 
single-qp current where the dynamical contribution of neutral modes is absent (cf. (HHj) ). 
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For weak renormalizatioii^ gc < /2v ((7^ < 5 for z/ = 2/5), as represented in figure [21(a), 
a maximum at ^ = V^ax separates the two above regimes, while for strong renormalization 
9c > p^/2z/ as in figure [2](b), the maximum is replaced by a smooth crossover a.tV = with 
a growing current. The two quantities V^^x defined at zero temperature as their 

notation remind us. In general both V^^x proportional to Un- However, their 

precise values depend on the parameters such as the filling factor and the renormalizations. 
For example, V^^^ can be obtained from fH6l) . in the limit Uc 00, by solving the equation 

iFi [2gn(3l 2ug,a\ + 2gn(3l - 1, -e*V^Jun] = . (51) 

For z/ = 2/5 and g-n = gc = 1, we find e*V^ax ~ 2.5 0;^ in accordance with the position of 
the maximum in figure [2]^a). 

At finite temperatures and low voltages e*V <^ T the current shows a linear behaviour 
~ G^^^V proportional to the linear conductance G^^\ The latter is represented in 
the insets of the figure [21 at temperatures smaller than the neutral bandwidth, T Un, 
neutral modes affect its scaling with G^^^ oc T'^('^^''°'^~^^"^i-~^\ On the other hand, at T ^ tu^ 
the scaling is driven by the charged mode only with G^^^ oc T'^('^^''°'^~^^ [32]. These two 
power-laws are distinguishable in the curves presented in the insets. 

At larger voltages the scaling behaviour fH9l) with neutral modes is still visible for voltages 
range T/e* -C <^ V^^x weak renormalizations (figure [21(a)) or for T/e* <^ -C 
for the strong case (figure [21(b)). Under these conditions the maximum and the 
crossover position remain substantially uneffected varying the temperature. Only at higher 
temperatures T/e* > V^^x T/e* > V^, depending on the renormalization, the neutral 
modes are saturated at any voltages. Their dynamical contributions disappear and the 
current directly pass from the linear behaviour to the power-law ( l50l) determined by the 
charged modes only. In these high temperatures condition, represented in figure [21 with 
orange and red lines, the position of the maximum Knax or the crossover Vc starts to increase 
linearly with temperature. 

The above discussion shows two important aspects: 

i) renormalization effects change qualitatively the current behaviour; 

ii) neutral modes dynamics may be visible in the current. 

In particular, we have seen that weak interactions may be most favourable in order to achieve 
point ii). Indeed, the experimental verification of a temperature independent maximum 
position, which would be a trade-mark of the presence of neutral modes is more easily 
detectable than the measurement of a temperature independent crossover that would be the 
case for strong renormalization. 

% Note that the classification of weak and strong renormahzations adopted here is hmited to the analysis 
of single-qp current. Later we will consider different definitions in order to follow the richer phenomenology 
introduced by the presence of the p-agglomerates. 
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We would like to conclude by commenting that if it seems not too difficult to give evidence 
of the neutral modes, much more delicate is the possibility to quantitative determine the 
neutral bandwidth since the current behaviour depends crucially on other a priori 
unknown parameters. 

4-2. Single-qps and p- agglomerates crossover phenomena 

In this part we will consider the behaviour of the total backscattering current ((/tot))iin 
( l43l) given by the current contributions of the two most relevant processes: the single-qp 
((j'-^-'))^, and the p-agglomerate ((J*-p^))^. The aim will be to investigate the crossover 
phenomena among these two different terms. We remind that the weight of these tunnelling 
processes are given by ti = t and \tp\ = Hp\t\. In the following Hp will be assumed as a free 
parameter. Since the experimental data, such as [30], shows the presence of single-qp at 
larger voltages/temperatures and the presence of p-agglomerates only for lower energies one 
need to choose Hp sufficiently low such that the effects of the p-agglomerate are not visible 
at high energies but it has to be big enough to see them at sufficiently low energies. 
In order to determine {{l^^^}}^ we remind that the most relevant p-agglomerate operator 
in ( TT5l) is given hj d = q = 0. This implies /3p = in iHM . with no presence of neutral 
modes in the dynamics. Only charged modes will enter with the coefficient = 1 as stated 
in ([9]). This allows to determine the current ((J*-^^))^ by using the results obtained in the 
previous section for the single-qp process with topological neutral modes 0;^ = - see (HTIl 
- apart from the necessary substitution ai ap = 1. 

Figures M^a) and [3](b) show the total backscattering current as a function of voltage for 
unrenormalized and renormalized cases respectively, and u = 2/5. 

We start the discussion at zero temperature corresponding to black lines. The current shows 
three different power-laws 

{y2^9.a?-i 1/ < 1/* (52) 

y2ug,ai+2gr.(3i-l < ^ < uj^/c* (53) 

V^'^^'''-^ LJn/e* < V, (54) 

with ctp, «! in and /?i in (1441) . Here, we defined with V* the bias at which the two main 
current contributions, at zero temperature, are equal ((-^^^^))^ = The precise 

value of this point crucially depends on Kp and on other parameters. 

Note that the power-laws (!52|) - (l5^ depend on the parameters gc^n which may affect 
qualitatively the behaviour. We already discuss their influence on the single-qp contribution, 
here we would like to focus mainly on the p-agglomerate part, which dominates the current 
at low voltages, see fl52l) . and on the crossover region. 

For weak values of charge renormalizations 1 < gc < 1/21/ the current decreases with voltages 
at V < V*, see (1521) and the black line in figure [3]^a) present a minumum around V*. This 
behaviour is completely different from the case of larger renormalizations gc > l/2z/, which 
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Figure 3. Total backscattering current {{Itot))i as a function of e*T^/aJoi with ujq = lO^^Uc, 
in bilogarithmic scale for = 2/5. (a) Unrenormalized case (7c = 5n = 1 with Kp — 0.4; (b) 
renormalized case gc = 5.5 and gn — 2, with Hp = 10. Temperatures are T/luq = (solid 
black), 0.025 (dot-dashed blue), 0.1 (dashed green), 2 (short dashed orange). The neutral 
mode bandwidth is u;„ = 10~^a;c. Unit of current: It = e|tp/(47r^a^Wc)- 



shows an increasing current (cf. figure Et^b)). In this last case, there is a tendency to 
lose the clear evidence of the p-agglomerate which merges very smoothly on the single-qp 
contribution. Indeed, the crossover is only signalled by a change in the power-law that is 
gradually less evident increasing charge renormalization. Hence, we can conclude that, at 
T = 0, the weak interaction case with 1 < gc < l/^^^, which exhibits a minimum in the 
current, is the most favourable in order to show the presence of agglomerates. 
We would like now to investigate the robustness of the above results for finite temperature. 
The main task is to extract the voltage and temperature ranges where one could clearly 
detect the contribution of p-agglomerates. At finite temperatures (colored lines in figure [3] 
), and low bias V -C T/pe* < V* the current, instead of the power- law (l52l) . shows a linear 
behaviour typical of the ohmic regime preventing divergence at very low voltages for weak 
interactions {gc < l/2z/). On the other hand at larger voltages T/pe* -C F < l^* it is 
possible to detect the presence of p-agglomerates from the behaviour of the current that 
follows the power-law (!52|) . At larger voltages V V*, single-qp dominates and current is 
similar to the T = cases (!53|) and (!54|) . 

As a common trend we observe that the optimal possibility to detect agglomerates depends 
crucially on the interaction parameters gc^n- From figure [3t^a), representative of weak 
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renormalization g^ < l/2z/, one can see that the transition between the hnear behaviour 
to the p-agglomerates scahng can create a peak around pe*V ~ T. In this case the current 
exhibits two peaks: the first due to the crossover between ohmic behaviour and the p- 
agglomerate power-law, the second due to the crossover inside the single-qp contribution 
due to the de-activation of neutral modes. Increasing interactions ( gc > l/2i^) the first 
peak tends to disappear until for strong enough interactions, gc > p^/2z/, as in figure [^b), 
even the second maximum disappears in accordance with the behaviour of the single-qp 
contribution discussed in the previous section. Here, the change in the power-law induced 
by the appearance/disappearance of p-agglomerates becomes less visible. 
The above discussions demonstrate that is possible to find voltages and temperatures regimes 
where the agglomerate contributions are detectable in the current. In order to confirm this 
important result we would like to find stronger signatures in other transport properties. The 
study of higher current moments, performed in the next section, will represent a powerful 
tool to clarify the crossover physics between the p-agglomerates and single-qps. 



4-3. Fano factor and higher cumulants 

Noise measurements are important tools in order to obtain clear information on the carrier 
charges involved in transport and consequently on the dominant excitations present in the 
tunnelling current. 

In order to study the current fluctuations we start by considering the Fano factor F2 = 
{{hot))2l^\ {{hot))i \ defined as the ratio between the zero frequency noise of the total 
backscattering current and the current itself. Using (138|) one can write the Fano in terms of 
single-qp and p-agglomerate current contributions 

e* ((j(^)))^coth(/je*V-/2) +p((jW))^coth(/3pe-l^/2) 

In the following we will restrict the discussion at energies smaller than the neutral 
bandwidth, i.e. e*V,T < Un-, for which crossover phenomena appear. Figure H] shows 
the Fano factor as a function of bias at z/ = 2/5 in the unrenormalized case, panel HJ^a), and 
in the strong renormalized case, panel |l](b). At zero temperature (black lines) the system 
is in the shot-noise regime, from (!55|) one then has 

_ e- ((/'")). +p ((/")). 

e ({/(■)», + {{/("», ■ ^""J 
We remind that the total current in flS2]) - flM|) is dominated by the p-agglomerates for V <^ 
V* where {{Itot))i ~ ((^^^^>>i, and by the single-qp at 1/ > where {{Itot))i ~ ((^^^^>>r 
Hence, the Fano shows two different plateaux 

^ = zy V ^V* (57) 

— = - y > K*. (58) 

e p 




10"'' 1 10^ 10"^ 1 10^ 

Figure 4. Fano factor F2 as a function of bias e*V/u)Q, with loq = 10~^Wc, in linear- 
logarithmic scale for v — 2/5 and for a;„ = 10~^Wc- (a) Unrenormalized case gc — g-n — 1, 
with Hp = 0.4; (b) renormalized case gc = 5.5 and gn = 2, with Kp — 10. Different line 
styles correspond to temperatures: T/luq = (solid black), 0.025 (dot-dashed blue), 0.1 
(dashed green), 0.5 (short dashed orange), 2 (dotted red). 



At low bias - cf. fl57|) - the Fano is dominated by p- agglomerates with effective charge pe*, 
while at higher voltages - cf. fl58l) - the single-qp prevails with charge e*. This behaviour is 
clearly showed in figures Hl^a) and IH^b) and is stable against renormalization effects. Low 
and high bias regimes are separated by an interpolating region. As shown in [32] the width 
of this region is not universal and crucially depends on the interaction strength. 
Note that the presence of two plateaux in the Fano is a direct evidence of the existence of 
two different kind of carriers that participate to tunnelling. At the zero temperature their 
charges are deduced by the values of the plateaux themselves. 

In order to compare with experiments we need to generalize the above results at finite 
temperatures. We will see that some previous conclusions are still valid despite the presence 
of more intricate regimes. At small voltages pe*V -C T the noise is represented by the 
Johnson-Nyquist thermal contribution with F2 = 2kBT/eV and is diverging for V ^ 0. 
This is clearly visible in all the finite temperature curves of figure HI For the discussion 
at larger voltages we consider separately weak and strong renormalizations. Let us start 
with figure 111(a). For temperatures smaller than the threshold voltage V* and voltages in 
the range T /pe* <V* the Fano reaches the first higher plateau aX F2 ^ v related to 

the charge of the p-agglomerate. On the other hand, at higher bias, it drops to the second 
plateau at F2 ~ u/p corresponding to the charge of the single-qp. 
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At larger temperatures T > e*V* the visibility of the p- agglomerate is compromised by the 
presence of thermal noise. Indeed, the Fano directly passes from the Johnson-Nyquist regime 
to the single-qp plateau (see e.g. orange and red curves in figure Hl^a)) without a possibility 
to detect the p-agglomerate plateau. One then can say that, for weak interaction one may 
detect the presence of two different excitations with a proper determination of charges only 
if T <C e*V*, that typically in the experiments correspond to very low temperatures. 
For strong renormalizations, as in figure Hl^b), the results are even more involved. While the 
extremely low and high temperature regimes with T e*V* and T ^ e*V* show a similar 
behaviour as observed above, at intermediate temperatures qualitatively differences appear. 
For example green line in figure Hl^b) shows two plateaux, but the first is at an intermediate 
value v/p < F2 < v smaller than the p-agglomerate charge. The origin of this behaviour 
is due to the smooth crossover between the p-agglomerate and the single-qp contributions 
present for strong interaction. Indeed, differently from the unrenormalized case, there is a 
finite region of bias where the two current contributions are comparable. Hence, the plateau 
value at Fg* corresponds to a weighted average of the single-qp (e*/e) and p-agglomerate 



We can then summarize the following important results: 

i) the existence of more than one plateau in the Fano factor is a clear trademark of the 
presence of more than one relevant tunnelling entity; 

ii) whenever more than one relevant tunnelling process is simultaneously present is not 
appropriate to evaluate the effective charge involved in tunnelling by fitting the noise 
in fl38l) using a single carrier charge as free parameter. 

This second procedure is often used in experiments [HI [151 ISSl [30], however it could bring 
to misinterpretations. For example, in the above case with an intermediate plateau at 
F2*, the current is still composed by two different excitations with different charges equally 
important. In this case fitting the Fano with a single effective charge will give a wrong 
description. All these observations could be relevant for the experiment in [30] especially 
for V = 3/7, where apparently the Fano factor instead to be equal to 3/7, even at low 
temperatures, reaches only the value of 2.4/7. 

In this last part of the section we would like to generalize the above discussion considering 
higher current cumulants. In view of their statistical independence, we expect to obtain 
complementary information on carrier charges from them. We define the generalized Fano 
factors as the normalized k-th order zero frequency current cumulants [53] 




with {{Itot))k the total cumulant defined in fH3l) . In the following we will consider the 
behaviour of odd cumulants with a detailed analysis of the normahzed skewness F3. This 
quantity is indeed experimentally accessible for a QPC in a Hall bar [ITj, [IS]. The general 



{pe*/e). 




{{hot))k 



(59) 
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Figure 5. Normalized skewness F3 as a function of bias e*V/ujQ, with ujq = lO^^cjc, in 
linear-logarithmic scale for i/ = 2/5 and cj„ = lO^^Wc- (a) Unrenormalized case gc = g„ = 1 
with Kp = 0.4; (b) renormalized case gc — 5.5 and gn = 2 with Kp = 10. The different 
line styles correspond to different temperatures according the same convention adopted in 
figure m The units are the same as figure [H 



expression for odd fc-th normalized cumulants can be easily obtained inserting fH2|) and (H3l) 
in (I59D 

e ) {{hot))i 
Note that this expression also holds at finite temperatures for odd k-th cumulant. For even k- 
th cumulants a coth(me*V/T) term appears such that ( 160|) is only valid at low temperatures 
T e*V. The absence of this term in the odd cumulants avoids divergences induced by the 
thermal regime which are, on the other hand, always present for even cumulants. Following 
similar analysis as done for the Fano we can extract the zero temperature behaviour of fl60l) 



.(0) 



pe 

e 

e* 
e 



k~l 



k-l 



. k-l 



k-l 



(61) 
(62) 



Again two plateaux appear with values related to the p-agglomerate charge at low voltages 
and to the single-qp charge at high voltages, with the power (A; — 1). The skewness in figure[5] 
reproduce this behaviour (black lines) both for the unrenormalized case of figure M^a) and 
for the strong renormalization of figure [S](b). 
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Increasing temperature the two plateaux are still present with a decrease of the higher with 
respect to the T = value. Note that the thermal divergence is not present here. Therefore, 
from the odd cumulants behaviour, one can have information on different excitations 
(looking at the different plateaux) even in the deep thermal regime e*V T which is 
not an observable regime in noise measurements. 

To be more specific let us compare the two red curves in figures [5t^a) and|l](a). While for the 
Fano the thermal regime cancels out completely the presence of the higher plateau, this one 
is still clearly visible in the skewness. This means that odd moments can be an important 
tool in order to detect the presence of different excitations (e.g. the p-agglomerate and the 
single-qp in our case) present in tunnelling processes. 

More delicate issue is the possibility to measure the precise values of the excitation charges. 
Indeed, as we saw before, one has to consider low enough temperatures (blue curves in 
figure [5]) in order to measure, from the higher plateau, the value of the charge associated 
with the p-agglomerate. This can be explained by observing that for — > 0, the current is 
proportional to the linear conductances (I4T1) . consequently, from (j60l) the cumulants are 

As extensively discussed in [32] one can verify that for temperatures low in comparison to the 
neutral modes bandwidth, T -C tu^, the temperature power- law of the single-qp conductance 
is C-^-* oc T^[^'='^/P^+'^^'*"^/p)^"~^] with positive exponent, while for the p-agglomerate one has 
G'^P^ oc T'^^^"'^"^^ with an exponent that changes sign for gc> l/u. 

Therefore for weak renormalization, g^ < l/i^, decreasing temperature, the p-agglomerate 
conductance strongly increase, and, at the same time, the single-qp contribution becomes 
negligible. Then one does not need extremely low temperatures in order to measure the 
true charge of p-agglomerate from the value of the plateau. This is well represented in 
figure [5](a) (green curve) where the skewness at low bias approaches the value v'^. Note 
that this charge would have been more difficult to be extracted from the Fano behaviour 
where the plateau is not so clearly defined - see the same curve in figure Hl^a) -. For strong 
renormalization > l/i/, as in figure [5](b), the power-law exponents of the conductances 
G'-^^ and G'^'^^ are both positives. Then, the dominance of the p-agglomerate appears at 
much lower temperatures with respect to the previous case and the plateau deviates from 
the T = value at temperatures lower than the weak renormalization case. In both case 
for temperatures not sufficiently low the intermediate values of the higher plateau depend 
on the weighted average of the two carrier contributions. 



5. Conclusion 



We investigated the effects of propagating neutral modes in the tunnelling of quasiparticles 
through a QPC for the Jain sequence of FQH. We observed that these modes qualitatively 
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modify the behaviour of the single-qp backscattering current. Even more interesting they 
determine the behaviour of the total backscattering current given by the sum of the single- 
qp and the p-agglomerate contribution. We demonstrated that the presence of neutral 
modes with finite bandwidth creates a crossover between the two kinds of quasiparticles 
involved in tunnelling. To avoid non universal effects that may hide the possibility to extract 
information on this crossover we considered the current fluctuations. These quantities are 
indeed more robust than the current in order to identify unequivocally the entities which 
tunnel. We first investigated the Fano factor that presents a plateaux structure related to the 
crossover between the single-qp and the p-agglomerate tunnelling regime. However, the Fano 
is strongly affected by thermal effects, therefore, we considered odd higher cumulants. We 
demonstrated that higher cumulants are an important tool in order to establish the number 
of relevant quasiparticle in the tunnelling. This is a strong motivation for investigating these 
quantities, despite of the enormous technical difficulties in their measurements. 
The cross-check of the physics described in this paper with some other independent 
transport quantities such as the finite frequency noise [Ml EHl [79] could open extremely 
interesting perspectives in the comprehension of anomalous behaviours that emerge in noise 
measurements. We are confident that analogous approaches could explain some unclear 
experimental observation for also the case u = 5/2. [STJ [80] . 
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Appendix A. 

This Appendix is devoted to the derivation of the coefficients Pm in ([141) . Let us start 
to discuss the monodromy condition introduced in section 12.11 The generalization of the 
internal statistical angle ([6]) for two different agglomerates with m and m' is the so called 
mutual statistical angle 9 defined as [61] 

¥"'\x)¥"''\x') = ¥"''\x')¥"'\x)e~' ® sgn(x-z')^ (A.l) 

the phase depends on the coefficients {a,P{k)) and {a',(3'{k')) that define two different 
operators in ([HD with (see ([9]) and ([TO]) ) 

a^ = j; pm{k) = + ^) + 2A;, (A.2) 

and similar expressions for am' and l3m'{k'). By using the operator commutation rules one 
obtains from flA.ip 

Q[am,Pm{k);am',Pm'ik')] = n [uamam' - Pmik)Prn'ik')] . (A. 3) 
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The monodromy condition imposes that the phase acquired by any excitations m in a loop 
around an electron must be a multiple of 27r [62],|63]. We then need, first of all, to start the 
discussion considering the electron operator, which corresponds to m = 2p + 1. For this we 
will generalize the results presented in [61], where the z/ = 2/5 case was discussed. In order 
to recover the electron operator of the original FL theory we consider f32p+i{k') in (]A.2p with 
k' = [61], we will check at the end the admissibility of this value. Other electron operators 



in (lA.ip with k ^ need to satisfy the mutual statistics constraint with the electron at 
k' = 0. From ([Ql) it is 

= TT [ua2p+ia2p+i - p2p+iik)P2p+iiO)] = Trr (A.4) 
with r G Z. Inserting flA.2l) in flA.4|) one has 

2 + -- (2 + -] J{p + l)(4p3 + 8p2 + 5j9 + 1 + 2kp) = r. (A.5) 
p \ pj 

To fulfil this equation one needs first of all to impose a constraint on the square root with 
{p + l)(4p^ + 8p^ + 5p+l + 2kp) = (A.6) 
Hence, the monodromy condition (1A.5P can be written as 

Ji^n-a)=r. A.7 
p 

One can see that to fulfil the two above equations one needs integer values of a. The more 
general solution of (1A.7I1 . is then written in terms of a number m G Z 

a = l + up, r = -u{2p + l). (A. 8) 

Replacing the above solution for a in (1A.6P one has 

pu"^ + 2u- (4p3 + 12p^ + 13p + 6) = 2k{p + 1). (A.9) 

This represents a parabolic Diophantine equation for u and k. In the following, we will 
restrict the analysis to the case with p < 6. Here, the only possible solution, is given by 

u = l + t{p + l), k=^{p+l){t-2){2 + 2p + pt), (A.IO) 

where t G Z. For p > 6 more complicated solutions are present and will be not discussed in 
this paper. Note that the t = 2 value in (lA.lOp gives /c = 0, in accordance with the initial 
choice A;' = in (lA.4p . Replacing the expression (lA.lOp for k in (lA.2p . with m = 2j9 + 1 we 
obtain 



«2p+l 



2 + -; P2p+i{t) = Vpip + T)(t + -]. (A.ll) 
p \ PJ 

Now we generalize the analysis to the m-agglomerates. The monodromy condition flA.SP 
between an m-agglomerate described by flA.2p and an electron, described by flA.lip is 

e _ m - ^y2kp + m'^{p+l)^y{p+l){pt + 1) _ 



TT p 



(A.12) 
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where f G Z. Similarly to the previous case, one has first to impose a constraint on the 
square root flA.12|) 

[2kp + m^{p + l)]{p+l) =b\ (A.13) 

The monodromy condition is then written as 

m — b(l + pt) . , X 

^——Jil = V. (A.14) 

p 

One can see that to fulfil the two above equations one first need integer values of b. The 
solution of (1A.14P is then 

b = m-fp, V = f + {fp~m)t ifeZ). (A.15) 
Replacing the above solution for b in (]A.13|) one has to solve the relation 

fp-2fm-m^{p + 2) = 2k{p+l), (A.16) 
whose solution can be easily derived, again for p < Q, 

f = -m- h{p + l), k = h{p+l){m+^p) (A.17) 

with k E Tj. 

Note that for m = 2p + 1 and choosing f = 2 — u one recover the equation flA.91) for the 
electron. Through the redefinitions m = sp + d and h = q — s {q E Z),we can rewrite ( ]A.17p 
as 

k = ^i^±ll (g2 _s^)+d{p+ 1) (q-s). (A.18) 

This corresponds to the expression f|T2l) quoted in the main text. Replacing flA.lSP in the 
initial expression ( 1A.2P for (3m{k) and using the notation m = sp + d one has 

= VpF+1) (q + ^) (A.19) 

as defined in (fT^ . Note that for s = 2, d = 1 and with the trivial identification q = t we 
recover the expressions for the electron in (lA.llll . 
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